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Abstract 

We study explicitly the applicability of perturbative QCD (pQCD) to the pion 
virtual Compton scattering. It is found that there are central-region singularities 
introduced by the QCD running coupling constant, in addition to the end-point 
singularities generally existed in other exclusive processes such as the pion form 
factor. We introduce a simple technique to evaluate the contributions from these 
singularities, so that we can arrive at a judgement that these contributions will 
be unharmful to the applicability of pQCD at certain energy scale, i.e., the "work 
point" which is defined to determine when pQCD is applicable to exclusive processes. 
The applicability of pQCD for different pion distribution amplitudes are explored 
in detail. We show that pQCD begins to work at 10 GeV 2 . If we relax our constraint 
to a weak sense, the work point may be as low as 4 GeV 2 . 
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The application of perturbative Quantum Chromodynamics (pQCD) to exclu- 
sive processes started since late 1970's [1,2,3], and now the pQCD approach 
has been widely employed to various exclusive processes. It is generally known 
that pQCD is applicable at very high energy scale, and one of the challenging 
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questions is whether pQCD is valid or not at present experimental accessible 
energy scale. A typical example is the applicability of pQCD to the pion form 
factor. It should be actually unknown whether pQCD is applicable or not be- 
fore any affirmative judgement that the dominance of perturbative effect over 
the soft effect is justified. It was indicated by Isgur and Llewellyn Smith [4] 
that a significant fraction of contributions to the form factor is from the soft 
end-point regions at medium-to-high energy scale. In calculating the pion elec- 
tromagnetic form factor, Huang and Shen [5] showed that the non-perturbative 
contribution from end-point regions can be suppressed with the end-point- 
suppressed distribution amplitude by taking into account the transverse mo- 
mentum in the pion wave function. Li and Sterman [6] reduced the end-point 
contribution by Sudakov effect, i.e., to replace the end-point singularity by 
theoretical based resummation over soft radiative corrections. There is also an 
attempt [7] to explain the large discrepancy between theory and experiment 
by large contribution from next-to-leading-order power corrections. 

One of the alternative exclusive processes, Compton scattering, which is an- 
other good laboratory to study the structure of hadrons, has been extended 
from real photon to virtual photon for the proton case as well as the pion case. 
Several independent theory groups have calculated the pion virtual Compton 
scattering (VCS) with very different methods, such as the chiral perturba- 
tion theory [8], the QCD sum rules [9], and pQCD approach [9,10,11,12]. One 
should meet the same question concerning the validity of pQCD, in similar to 
the case of the pion form factor. It is the purpose of this paper to study explic- 
itly the applicability of pQCD to the pion VCS process. The early literature 
focused on how to eliminate the specific problem of kinematic singularities by 
analytic integration, but lack of a full study of the end-point effect. We will 
show that there are central- region singularities in the pion VCS, in additional 
to the end-point singularities generally existed in other exclusive processes 
such as the pion form factor. We will introduce a technique to evaluate the 
contributions from these singularities, so that we can have a criterion for the 
applicability of pQCD at certain energy scale. 

The pion VCS process is illustrated in Fig. 1, where and are the 4- 
momenta of the incoming pion and virtual photon respectively, and p" 1 and 
q'^ are the corresponding 4-momenta of the outgoing pion and photon. We 
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define the following frame invariant variables: 



Q 2 = -q 2 , S=(q + p) 2 = (q' + p') 2 , 



(1) 



t= (p'-p) 2 , 



u=(p' - q) 2 , 



where q 2 is the squared 4-momentum of the virtual photon, and S is the 
squared center-of-mass energy of the virtual photon and pion system. 



Fig. 1. Virtual Compton scattering 7T7* — > 7T7 shown in the center-of-mass frame. 

Factorization theorem [1,2] tells us that the scattering amplitude of an exclu- 
sive process can be expressed in a convolution formula: 



M(Q 2 )= J[dx] J[dy}cl ) *(y h Q y )H(x h y i ,Q 2 ) ( l ) (x i: Q x )[l + 0(rn 2 /Q%(2) 



where in the pion VCS case [dx] = dxdx5(l — x — x), x = 1 — x, Q x = 
mm(xQ,xQ), x is the light-cone k + = k° + k 3 momentum fraction of the 
struck quark (antiquark) in the incoming pion, and y is the light-cone mo- 
mentum fraction of the corresponding quark (antiquark) in the outgoing pion. 
The validity of Eq. (2) is that there is an initial Q = Q as the factoriza- 
tion scale above which the hard part and the soft part could be separated 
safely and explicitly. We expect S, —t, and — u large to guarantee factoriza- 
tion legal in the VCS case. The distribution amplitude <j>(x, Q) may be chosen 
without higher Gegenbauer polynomials to exclude Q dependence. In this pa- 
per the asymptotic (as) [1,2] and Chernyak-Zhitnitsky (cz) [13] distribution 
amplitudes: 




i 



i 



o 



o 
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<f> as = y/3f ir x(l-x), (3) 
<p cz = 5V3fnx(l-x)(l -2xf, (4) 

and their end-point-suppressed partners of Brodsky-Huang-Lepage (bhl) [3,14] 
and Huang-Shen (hs) [5] distribution amplitudes: 



<t> bhl = 1.4706^/^(1 - x) exp [-^S, (5) 

X I X X J 

<j) hs = 8.8763^/^(1 - x)(l - 2xf exp f"' ™ 62 ], (6) 

X I X X ) 

are used and compared in our calculations, where the pion decay constant 
f n = 93 MeV is adopted. We use the leading-twist factorization scheme instead 
of the handbag scheme. The different hard scattering amplitudes H have the 
following relations due to parity symmetry [12]: 



Hll — H RR , Hrl — Hlr, 

H+l = H +R , H- R = —v^H+h, H L = —v^H+l, 



thus only 3 independent H LR , H RR , H +R needed to be evaluated. The expres- 
sions of the hard scattering amplitudes H' [12] , which have been multiplied 
by the prefactor xxyy of distribution amplitudes and not included ^a e a s yet, 
are listed in Table 1, where the notations for the diagrams are corresponding 
to Fig. 2 of Ref. [12]. 

For simplicity, we neglect the mass of pion and work in the center-of-mass 
frame (Fig. 1) following available pQCD studies [9,10,11,12]. The 4-momenta 
of the incoming photon and pion are: 

q» = (u,0,0,p), p"=(p,0,0,-p), 
and the corresponding 4-momenta of the outgoing particles are: 



g ,M = ^±^(1, sin0,O, cos 0), p ,M = ^4^(1, - sin 6, 0, - cos0). 
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Table 1 

The hard scattering amplitudes calculated by pQCD 



diagram 



LR 



H't 



RR 



+R 



20xy 
9v(x—a) 



29xyc 2 
9vs 2 (x—a) 



AOcxy 
9s 



20x 
9v 



20xcf_ 

9vs 2 



20cx 
9vs 



-2Qxc z 



9vs 2 (x—a) 



-20cx 
9s 



16xx ys 2 i 1 1 n 

9c? \x—a x—b> 



16xy(l-vxs 2 ) ( 1 1 



9vc 2 



\x— a x—b> 



— 16sxyy(l—2vxs 2 ) 
9vc(l—ys' 2 )(x—b) 



-Wxxs 2 I 1 _ 1 \ 
9c? ^x— a x—b' 



IGcsxxy 
9(l-ys 2 )(x-b) 



20y 

9vs 2 



20y(l-2vx) 
9vcs 



—20axy 
9[y(l-vs 2 )-v] 



-20xyc 2 
9vs 2 [y(l — vs 2 ) — v] 



—20caxy 
9s[y(l— vs 2 ) — v] 



-20y[y+vx{l-2ys 2 )} 
9v[y(l—vs 2 )—v] 



~20y[l-s 2 (y+vx)+2vxys 4 )] 
9vs 2 [y(l — vs 2 )— v] 



-20cy 2 (l-2yxs 2 ) 
9vs[y(l— vs 2 )— v] 



— 16xxc 2 s 2 y 2 

9(l-ys 2 )[y(l-vs 2 )-v](x-b) 



— 16xxyc 2 

9[y(l— vs 2 )— v](x- 



— Wsxxy 2 c 3 

9(1— ys 2 )[y(l— vs 2 )— v](x- 



— 16xyys 2 [y— vx(l— 2ys 2 )] — 16xyys 2 [y— v+vx(l— 2ys 2 )] — 16csxyy[v— y+2vxys 2 } 

9v(l—ys 2 )\y(l—vs 2 )—v](x—b) 9v(l— ys 2 )\y(l— vs 2 )— v](x— b) 9v(l—ys 2 )\y(l—vs 2 )—v](x—b) 



Here we define 



S = (<f + ^) 2 = (uj + p) 2 , 



c = cos ^, s = sin ij, 



v = q 2 /S, v = 1 — v, 



(8) 



a — _H l) — y-v-yvs 2 

v ' v(l— ys 2 ) ' 

t = {p'^-p") 2 = -Svs 2 , 
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where v stands for the photon virtuality and should be between —1 and 0, and 
6 is the angle between the incoming virtual photon and the outgoing photon 
and can be obtained from the experimental variables by 



We let v = —0.8 as an example in our calculations to ensure q 2 large, —t is 
monotonously increasing with the scattering angle 9 if S and v are both fixed. 
For diagrams a, b, c, f, g and h (see Fig. 2 of [12]), the squared 4-momentum 
transfer of the exchanged gluons between the two valence partons in the pion 
is 

Q 2 = -vxyS, (10) 



where the small Q 2 region is also the end-point region whose effects will be 
highly suppressed by the distribution amplitudes. For diagrams d, e, i and j, 
where the incident and outgoing photons are connected to different quark lines, 
the corresponding squared 4-momentum transfer of the exchanged gluons is 

Q 2 = v(l-ys 2 )(x-b)S. (11) 



Because of0<a<6<l, whatever the kinematical region we will choose, it 
is impossible to guarantee Q 2 large enough for making pQCD legal here. The 
additional singularity of Q 2 in the central region, i.e. Q 2 — > when x — > b, 
makes calculations here much more complex than that of form factor. We use 
principle integration formula as in [12] 



l im f IM dx = P [ l^Ldx-mf(a) 
e->o J x — a + ie J x — a 

o o 

} r f( x ) — f(o-) p/ \ pi I — a . , ,„^ N 

= /{ — V - + / (a log m]}dx (12) 

J x — a a 
o 

to get reliable numerical results. 
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Straightly, we have: 

da 1 ^ tJfl2 3(2nfal(V3f n : 



2 



dcos9 32nS(l-v*) mi ' 2 -i K "''s^l-vi) J mf > {U) 

where £|M'| 2 = \M' LR + M' RR + (1 - l/v)M' +R \ 2 , with M"s calculated from 
Eq. (2) with if"s in Table 1 by including the factor a s (Q 2 ) which is not chosen 
as a frozen coupling constant [12] in this paper. M RR contributes dominantly 
over M\ R and M 2 R . For the applicability of pQCD, the concerned gluon should 
be hard, i.e., Q 2 should be large enough to prevent mis-absorbing soft contri- 
butions. Therefore, we do not fix the QCD running coupling constant 

a ^ 2) = h 1 %/A2 ' (14) 

b log Q 2 /A 2 

where bo = 11 — 2n//3, n/ is the number of active quark flavors, and A ~ 
200 MeV. We let a s running if Q 2 is larger than a given g GeV 2 ; otherwise, 
we let a s = 0, expecting that only pure pQCD contributions are considered. 
a s (g) and factorization scale Q can not be determined a priori by pQCD 
theory. So the choosing of g seems of some kind of arbitrary. However, the 
physics should not depend on which cutoff we make, i.e., the results should 
not change much if we change g slightly. In Fig. 2 we find that in the <f) as 
case different g do not affect the result through 10° — » 80° and deviate a little 
in 80° — > 170°. In the <pbhi case the two curves of different g are even more 
adjacent because <pbhi is just the end-point-suppressed partner of <f) as . <p cz has 
bumps near the end points so its results are very sensitive to g (Fig. 3). It 
can be expected that <phs will be better than <p cz at the cost of suppressing 
the end-points, just as we show in Fig. 2. With the squared center-of-mass 
energy S increasing, the gluons have generally more tendency to be hard. 
Thus, the parameter g becomes less important. But <p cz is still bad even when 
S reaches 50 GeV 2 , therefore we will not consider it any more as it lacks 
predictive power. We calculate the cross sections by using different distribution 
amplitudes and varying g = 0.09 (a max = 1.7), g = 0.49 (a max = 0.56) and 
S = 4,10,20,50 GeV 2 (Figs. 2, 4, 5, and 6). Fortunately, the cross sections 
are independent to g when the energy scale is large enough. This leads us to 
put forward "work point" , a concept analogical to "Kondo temperature" . 

When the cross sections of one given distribution amplitude are no longer 
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Fig. 2. Cross sections at S = 4 GeV 2 . 
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Fig. 3. Cross sections of (f) c 
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Fig. 4. Cross sections at S = 10 GeV 2 . 

affected by g, we define the present energy scale as the "work point" of this 
distribution amplitude. Different distribution amplitudes have different work 
points. The end-point-suppressed one usually has lower work point. Only above 
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Fig. 5. Cross sections at S = 20 GeV 2 . 
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Fig. 6. Cross sections at S = 50 GeV 2 . 

the work point which means that the contribution from Q 2 < g GeV 2 region 
can be safely neglected and all the contribution comes mainly from Q 2 large 
region, the condition for the applicability of pQCD is satisfied and the cal- 
culations of the given process using this distribution amplitude are possibly 
reliable. By our work it is clear that the work points defined by S are, as : 
~ 20 GeV 2 ; <p bM : ~ 10 GeV 2 ; <p hs : ~ 20 GeV 2 . With the transverse momentum 
contributions considered, Sudakov effect [9] may reduce the work point lower. 
If we relax our constraint to a weak sense, i.e., demanding the cross section 
difference with different g should be smaller than some percents, the work 
point of (fibhi and (p as will be lowered to 4 GeV 2 in the small angle region for 
Fig. 2. This is in agreement with Coriano-Li [9]. 

As shown in Fig. 7, we calculate S 3 d ^ s9 with S varying from 4 GeV 2 to 
50 GeV 2 . it^-b obeys approximately the S~ 3 law at fixed angle as naive pQCD 
has predicted [15]. The cutoff of the running a s suppresses the value of the 
cross section with approximate one order of magnitude, but it renders enough 
to survive the scaling law. It is also nontrivial that the cross sections are almost 
the same at small angles, especially from 10° to 80°, for different distribution 
amplitudes (Fig. 8), as we have no justification to demand that. On the other 
hand, they are so apart at large angles, so that with further experiments we 
can determine which one should be chosen in the VCS case. We should pay 
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Table 2 

Total cross sections [20°- 80°] 



S(GeV 2 ) 4 10 20 50 



a (nb):(f) bh i 


377 16.9 


1.75 


9.22 x 10~ 2 


a (nb):<^ s 


6.5 


1.52 


1.12 x 1Q- 1 



attention that the cross sections of (phs have a minimum running from 60° to 
70° while the cross sections of <j) as and <p bh i have the same minimum always 
at 80°. Summing up 20° to 80° to get the "total" cross section should be 
meaningful as predictions, as shown in Table 2. 

We should state that this is not an exact numerical estimate of the cross sec- 
tions, but rather a consistency check, without any higher-twist contributions 
or radiative corrections. SELEX Collaboration at Fermilab [16] obtained the 
total forward cross section of 7r~e — > n^e'-f of 38.8 nb, in agreement with the 
theoretical expectation by the chiral perturbation theory method [8]. Their 
S < M 2 , which is below the work point of any distribution amplitude we have 
explored. It illuminates that in the nowadays experimental available energy 
region the non-perturbative contribution is dominant over the pQCD contri- 
bution for the pion VCS process. 

The dependence on the photon virtuality v is not very sensitive, as can be seen 
from Fig. 9, where two different values of v = —0.3 and —0.5 are examined for 
(fibhi at S = 4 GeV 2 . In other words, the squared virtual photon 4- momentum 
transfer q 2 = Sv is not a signal energy scale in virtual Compton scattering. 
And S = 1 and 2 GeV 2 for (p^hi are also examined in Fig. 10, from which 
one may even arrive at a conclusion for the applicability of pQCD, in a rather 
weak sense, for 0^ at S — 1 — > 2 GeV 2 in the region 20° to 80°. 

Every distribution amplitude evolves in principle to the asymptotic form at 
very high energy scale. Such evolution effects are small for (p as and (pbhi, 
whereas they should be significant for <p cz and 4>hs ■ However, there has been 
no numerical evaluation on the evolution of (p cz and <p hs in previous studies. 
It is reasonable to expect a lower energy scale for the applicability of pQCD 
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Fig. 7. Cross sections at different S - scaling (g = 0.49). 
for cf) cz and <ph s than what we claimed, after taking into the evolution effects. 

In summary, we studied explicitly, for the first time, the applicability of pQCD 
to the pion virtual Compton scattering. As in other exclusive processes such 
as the pion form factor, the end-point contribution is not negligible in the 
pQCD calculation of the pion VCS process, as manifested by the difference 
between the cross sections of <p as with ip^i and <ft cz with (f>hs- We also noticed 
that there exist middle-region singularities introduced by the QCD running 
coupling constant of the exchanged gluons between the two valence partons 
of the pion. We introduced a simple technique to evaluate the contributions 
from these singularities, so that we can arrive at a judgement that these con- 
tributions will be unharmful to the applicability of pQCD at certain energy 
scale, i.e., the work point to guarantee the safety of pQCD. The work points 
for different distribution amplitudes are explored in detail in this work. For 
the end-point-suppressed distribution amplitudes such as (fibhi, we showed that 
pQCD begins to work at 10 GeV 2 . If we relax our constraint to a weak sense, 
the work point may be as low as 4 GeV 2 . 
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Fig. 8. Comparison between different distribution amplitudes (g = 0.49). 
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Fig. 9. Cross sections with different v. 
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